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TECHNICAL NOTE 3278 

ATTENUATION IN A SHOCK TUBE DUE TO UNSTEADY-BOUNDARY-LAYER ACTION 

By Harold Mirels 


SUMMARY 

A method is presented for obtaining the attenuation of a shock wave 
in a shock tube due to the unsteady boundary layer along the shock-tube 
walls. It is assumed that the boundary layer is thin relative to the 
tube diameter and induces one-d im ensional longitudinal pressure waves 
whose strength is proportional to the vertical velocity at the edge of 
the boundary layer. The contributions of the various regions in a shock 
tube to shock attenuation are indicated. 

The method is shown to be in reasonably good agreement with exist- 
ing experimental data. 


INTRODUCTION 

A shock tube consists of a fluid at high pressure (region 4 of fig. 
1(a)) separated by a diaphragm, from a fluid at low pressure (region l) . 
When the diaphragm bursts, a shock wave propagates into region 1 while 
an expansion wave propagates into region 4. A time-distance plot of 
these waves under ideal conditions is indicated in figure' l(b ) . Regions 
2 and 3 have the same velocity and pressure but have different tempera- 
tures. The interface between regions 2 and 3 is referred to as a con- 
tact surface. The analysis of the flow for perfect fluids is straight- 
forward (see, for example, ref. l). In an actual shock tube, however, 
viscosity and heat conduction can not be ignored. These lead to a bound- 
ary layer along the walls of the shock tube as indicated in figure l(c). 
The boundary layer introduces nonuniformities into the shock tube. Ana- 
lytical studies of this boundary layer are presented in references 2 to 
6. One of the important consequences of the wall boundary layer is that 
it generates weak pressure waves which catch up with and attenuate the 
shock wave propagating into region 1. This attenuation has been studied 
experimentally and analytically in the work of references 1, 4, 5, and 
6, and is the subject of the present report. It is assumed that the 
boundary layer is thin relative to the shock-tube diameter. This is a 
practical restriction, since most shock tubes are designed so that the 
core of potential flow is relatively uniform in order to permit aero- 
dynamic tests . 
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A few remarks concerning previous shock-wave -attenuation analyses 
are appropriate. In reference 4, the coordinate system is defined so 
that the shock wave is stationary. The flow between the shock wave and 
the contact surface is considered as a one -d imen sional steady flow so 
that at each instant the mass flow through the shock wave equals the 
mass flow at the contact surface. If the mass flow at the contact sur- 
face is known at each instant, the corresponding shock strength can he 
found. The mass flow at the contact surface is determined from the lo- 
cal boundary-layer displacement thickness and free-stream conditions 
corresponding to an unattenuated shock. However, it can he shown that 
free-stream conditions do not remain constant at the contact surface 
(because of perturbations induced by the boundary layer) . Moreover, the 
method of reference 4 does not take into account the existence of weak 
pressure waves between the contact surface and the shock wave. Since it 
is precisely these pressure waves which are responsible for shock atten- 
uation, the method of reference 4 can not be expected to give accurate 
quantitative results. It does, however, indicate some of the important 
parameters involved. Reference 1 uses the method of reference 4 in its 
studies of shock attenuation. At best', qualitative agreement with ex- 
periment is indicated. 

In reference 6, the flow perturbations due to the boundary layer 
are considered in detail. At each section of the shock tube, the veloc- 
ity and temperature variations associated with the boundary layer are 
averaged across the tube to provide an equivalent one -dimens ional flow. 
The wall shear and heat addition due to dissipation and heat transfer 
at the walls (all found from boundary-layer theory) are assumed to act 
on this equivalent one-dimensional flow, their action resulting in the 
generation of one -dimens ional pressure waves propagating in both the 
upstream and downstream directions. By integrating along characteris- 
tic lines, the attenuation of the shock propagating into region 1 is 
then found. The theoretical trends appeared to be in good agreement 
with the particular experimental results reported therein. 

Reference 5 considers the flow in a coordinate system which is sta- 
tionary with respect to the shock. The unsteady nature of the flow be- 
tween the shock wave and the contact discontinuity is associated with 
the receding of the contact surface with respect to the shock wave. 

Weak pressure waves are assumed to be generated by the recession of the 
contact surface. The magnitude of these waves is obtained by a one- 
dimensional averaging procedure similar to that of reference 6. These 
waves overtake the shock and result in attenuation. The method of ref- 
erence 5 ignores the contribution of region 3 and can not be expected 
to yield good quantitative agreement with experiment. 

Of the previous reports on shock attenuation, reference 6 appears 
to give the best agre ement with experiment . However, it can be shown 
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that the method of reference 6 does not give a completely valid repre- 
sentation of the wave phenomena induced by the boundary layer along the 
shock-tube wall. The deficiency of reference 6 is mainly associated 
with its use of wall shear in the determination of the perturbation 
pressure waves generated by the wall boundary layer. It is well recog- 
nized that. the effect of a boundary layer on its external flow is di- 
rectly related to the vertical velocity at the edge of the boundary 
layer. For example, reference 7, which is concerned with the Rayleigh 
(impulsive -plate) problem for a compressible fluid, shows that the bound- 
ary layer generates pressure waves in the externa], flow which are equi- 
valent to those which are produced if the wa ll moves normal to itself 
with a velocity equal to the vertical velocity at the edge of the bound- 
ary layer. The nonzero pressure gradient over a flat plate moving at 
high speeds (because of the finite displacement thickness of the boundary 
layer) is an equivalent steady -flow phenomenon (e.g., ref. 8). Thus a 
proper my to find the waves generated by the unsteady wall boundary 
layer in a shock tube is to base the calculation on the vertical velocity 
at the edge of the boundary layer. Such an analysis is presented herein. 
The quantitative results for shock attenuation thus obtained would be 
expected to differ from those of reference 6. The comparison between 
the method of reference 6 and that of the present report is discussed 
more fully in the main body of the report and in the appendix titled 
GENERATION OF PRESSURE WAVES BY WALL SHEAR AND HEAT ADDITION. 


ANALYSIS 

One-dimensional flow with mass sources is treated, and the waves 
generated by these sources are derived. The waves generated by unsteady- 
boundary-layer action in a tube are then found, assuming the process can 
be considered as a one-dimensional unsteady flow. The application to 
the shock-tube problem is then indicated. 


Generation of Waves by Mass Sources 


Consider one -dimensional uniform flow in a tube of constant cross- 
sectional area. Assume that weak mass sources, uniformly distributed 
across each cross section, are present and perturb the flow. Denote the 
net perturbation of a quantity (from the uniform-flow conditions) by A. 
(See appendix A for def inition of symbols . ) The equations of motion are 


p (lr + u ir) “ * (Momentum )^ 

l^ + p lr + u lr = m ( Conttouit y) > 


CD 


V 

Zp = a^Ap (isentropy) 
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where m = m(x,t) is the rate of mass addition per unit cross-sectional 
area per unit x. The presence of the sources generates waves. Let the 
superscripts + and - indicate perturbations associated with waves 
moving in the +x- and -x-directions, respectively. If £/c designates 
integration variables for x,t, the solution for the net perturbation at 
any point x,t may be expressed as 

Ap = Ap + + Ap 

Au = Au + + Au" 

= ^ (^> + - AP~) 

Ap = Zip 
bT 

where 

* 2T1TMT jf “(«>* - iv£) ^ 

4-00 

As " = 2l3VTtJ^£ 4'*'^) 45 

The integrations are conducted along the characteristic lines 
t = t - £x-£)/(a+u) and T = t - (£-x)/(a-u) in the plane. The 

upper limit on the integral for Ap" Is +® or -® depending on 
whether M< 1 or M > 1, respectively. Equation (2b) incorporates 
the acoustic relations Ap + = pa Au + and Ap - = -pa Au“. 


(2b) 

(2c) 


(2a) 


Generation of Waves by Unsteady -Boundary -Layer Action 

Consider a tube of uniform, cross section to have flexible walls 
such that a small normal velocity v can be generated at the walls. 

This is equivalent to mass entering the tube at the rate pv dZ, per 

unit x, where the integral is taken around the perimeter of the tube 
cross-sectional area. If the flow in the tube is considered as one 
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dimensional,^" the equivalent source strength is 

M-fiVJ 

4pv 
d 

•where d= 4A/7 is the' hydraulic diameter, and v = v(x,t) . The ex- 
pression for Ap can then he written as 

f - S [rTMj^ ^ - rrf) ^ + rru £ ^ 

00 

Similar expressions can he written for Au and Ap. In the case of 
waves induced hy houndary- layer action, the v in equation (4) refers 
to the normal velocity at the edge of the boundary layer. Note that a 
positive v results in compression waves, while a negative v results 
in expansion waves. 


^If the flow through the tube cannot he considered as one dimen- 
sional, it is necessary to consider each element of tube surface as an 
elemental wave source of strength proportional to the local value of v. 
The net wave strength at any point in the tube is found from an integra- 
tion over the entire tube surface. 

In the problem of attenuation in a shock tube, the shock wave is 
considered to be uniform, laterally, at each value of x, but to de- 
crease in strength with increases in x. Hence, a one-dimensional ana- 
lysis is permissable. Actually, small lateral variations of shock 
strength exist, particularly near the walls, since the boundary layer 
immediately behind the shock exerts a three-dimensional effect on shock 
strength. These lateral variations are ignored herein. 

The shape of the shock wave was previously studied (ref. 9). It 
was assumed that region 2 was infinite in extent (i.e., region 3 was 
neglected), and the problem was considered as a steady-flow problem in 
a coordinate system moving with the speed of the undisturbed shock. It 
was found that the shock wave assumes a parabolic shape, but that appre- 
ciable shock curvature is restricted to a region near the wall less than 
a boundary -layer thickness in extent. 




6 


NACA TN 3278 


ATTENUATION IN A SHOCK TUBE 

Equation (4) can be applied to find the attenuation in a shock 
tube. The details of the analysis are described in appendixes B to F. 
The resulting formulas for shock attenuation are presented in the fol- 
lowing section. The limiting solution for weak shocks and a numerical 
solution for an air-air shock tube are also indicated. 


Shock Attenuation Formulas 


The flow in a shock tube is assumed to consist of the ideal basic 
flow plus small perturbations due to the boundary layer. Ideal shock- 
tube flow relations are summarized in appendix G. The expansion wave 
of the ideal flow is assumed to have negligible thickness and to propa- 
gate into region 4 with velocity u = - a^ as indicated in figure 2. 
Let point d of figure 2 represent an arbitrary point on the shock-wave 
characteristic. The problem is to find the net pressure perturbation 
behind the shock (i.e., Ap 2 This requires an integration of equa- 

tion (4) along a l 1 the characteristic lines which contribute to Apg^. 



be, and ab. Hence, these are the characteristic lines considered in the 
present analysis (as is indicated in appendix B) . Integrating along 
these characteristic lines permits dpg ^ to be expressed as 


1 a 2 d A P2,d 

F 2y £ p 2 




( 5 ) 


Equation (5) is derived in appendix F. The notation is defined in 
appendixes A and F. The term CE Apg c /4p 2 a represents the contribu- 
tion to Ap 2 ^ of all the characteristic lines other than bd, be, and 
ab (appendix F ) . 


If the boundary layer in region 2 is wholly laminar or wholly tur- 
bulent, the boundary -layer solution for v 2 is known (appendixes D and 
E). Similarly, if the boundary layer in region 3 is wholly laminar or 
wholly turbulent, the solution for v 3 is known (appendixes D and E). 
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For these special cases, the integrals of equation (5) can he readily- 
evaluated. Thus, if the houndary layers in regions 2 and 3 are hoth 
laminar, equation (5) becomes (from eqs. (F6) and (F7) with 

**2 = n 3 = V 2 ) 




' a 2 d ^2,d 

2 p 2 


-irgLgF 


fpT 5 


1 + 


( u s/ a 2)' 


Mg 




^52^2 

a 32 



C u s / a 2 ) 

M 2 


1 + M3 + a 43 



( 6 ) 


where Lg and L 3 
in regions 2 and 3 
(F6) and (F7) with 




are found from appendix D. If the boundary layers 
are hoth turbulent, equation (5) becomes (from eqs. 


n 2 


= n. 


= 1/5) 


(?)' 


—2 . 50 “f 1 

“s \ 
^ a 2/ 

i 

n . 1 - <u s 

/ a 2) 

1/5 

IT “2 . 



1 + 1 ~ ^ u s/ a 2 ^ 

p T52 M 2 L Mg 

a 32 1 + M 3 + a 43 



( 7 ) 


where Lg and L 3 are found from appendix E. Similar expressions can 

be found for the mixed cases ng = l/2, n 3 = l/5, and ng = l/5, 

n 3 = l/2. For other boundary-layer characteristics (such as transition 

from laminar to turbulent flow in the middle of regions 2 or 3), it is 
necessary to integrate equation (5) with the corresponding v distri- 
bution. Reliable criteria for determining the transition points in 
shock-tube boundary layers have not yet been established. A crude ten- 
tative method for estimating the transition points in regions 2 and 3 
is presented in appendix H. It is pointed out therein that the Reynolds 
numbers at point b, computed separately for regions 2 and 3, can be -used 
as an index to dete rmin e whether the boundary layers in regions 2 and 3 
are primarily laminar or turbulent, respectively. The values of Re 2 ,b 

and Re 3 as defined in appendix H, can be obtained from figure 3 for 
the case where the gas in region 1 is air at 520° R. 
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The present theory requires that the boundary layer he thin rela- 
tive to the shock -tube diameter. The houndary-layer thickness in region 
2 is defined hy equations (D5) and (E3) for 1 aminar and turbulent bound- 
ary layers, respectively. The laminar houndary-layer thickness 8g is 

taken to correspond to u/ug = 0.99, •while the turbulent boundary -layer . 
thickness Sg is taken to be the value obtained from, a Karman-Pohlhausen- 
type integral solution. Considering conditions along characteristic 
lines, the maximum, boundary-layer thickness occurs at point b of region 
2. Values of 5g and §2 ; b are plotted as a function of Mg in 
figure 4 for the case -where the gas in region 1 is air at 520° R. 

It is previously noted, in connection with equation (i), that a 
positive v results in the generation of pressure waves, while a nega- 
tive v results in the generation of expansion waves. From the 
boundary-layer theory of appendixes D and E it can be seen that, apart 
from dissipation and heat-transfer effects, Vg is negative and V3 is 
positive. Thus, the boundary layer in region 2 induces expansion waves 
(which attenuate the shock), while the boundary layer in region 3 induces 
compression waves (which accelerate the shock) . Dissipation and heat 
transfer modify these results. Dissipation tends to increase v in 
both regions 2 and 3. Heat transfer from the wall to the boundary-layer 
increases v. Heat transfer from the boundary layer to the wall de- 
creases- v. In region 2 the heat transfer is from the boundary layer 
to the mil (ref. 2) and leads to larger negative values of vg and 

therefore more shock attenuation. In region 3 the heat transfer is from 
the wall to the fluid for a weak expansion wave and from the fluid to 
the wall for a strong expansion wave (ref. 2). Thus, for a weak expan- 
sion wave, heat transfer in region 3 tends to generate compression waves 
(thus accelerating the shock) . For a strong expansion wave, the heat 
transfer in region 3 tends to induce expansion waves (which attenuate 
the shock) . The net effect of all these factors is to attenuate the 
shock wave . The relative magnitudes of the various terms in equations 
(6) and (7) are noted later in a numerical example. 

Limiting Solution for Weak Shocks 

For Mg approaching 1, equations (6) and (7) take an the follow- 
ing forms: 


C for 

(8a) 

( for T 2 ,w = T 2,r> 

(8b) 



ficn=& 
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Turbulent case: 

C for T 2,w = T l> 

(9a) 

(for T^ w = T 2 ^ r ) 
(9b) 

The condition T 2 v = Tj (eqs. (8a) and (9a)) corresponds to the case 

where the shock-tube wall is a perfect conductor (assuming that the 
fluid in region 1 is in thermal equilibrium with the wall) . Note that 
taking T 2 equal to Tp means that the wall remains at its original 

temperature and that heat is transferred from the fluid (of region 2) 
to the wall. The condition T 2 w = T 2 r ( e< l s * (db) a nd (9b)) corres- 
ponds to the case wherein the wall is a perfect insulator or has a very 
low heat capacity. When the shock-tube wall is a metal, the assumption 
T 2 ¥ = T : should give very accurate results, as discussed in reference 

2 (particularly since weak shocks are now being considered) . Results 
for both T 2 w = Tp and T 2 w = T 2 r are g iven 80 as to define the 
upper and lower bounds of the heat -transfer effect on shock attenuation. 
The attenuation is greater, by a factor of approxiunately y 1 , when there 

is heat transfer as compared with the insulated-wall case. (The factor 
is exactly Tp for the turbulent case and approximately Tq for the 

laminar case, the latter depending on the value of ffq) . 


pj / 5 ^agdj -/ 5 Ap 2f 


*2 


-o.nsrp 
[zCrp + i)] 
-o.iisrp 
[z(ri + i)] 


1/5 


1/5 


Air -Air Shock Tube with T^ = Tp = 520° R 

Equations (6) and (7) were evaluated numerically for an air -air 
shock tube with T^ = Tp = 520° R. It was assumed that a = 0.70 and 

y = 1.4. The results for the laminar and turbulent cases are summarized 
in figure 5. Both the insulated-wall case (T g = Tg r J T 3 w = T^ r ) 

and the case where the wall is a perfect conductor (T 2 = T 3 v = Tp =? T^_ 

are noted therein. When the shock-tube wall is a metal, the assumption 
that the wall is a perfect conductor should give accurate results ex- 
cept possibly for very strong waves (e.g., ref. 2). 


The relative contribution to the net attenuation of the various 

terms in equations (6) and ( 7 ) is indicated in figure 6 for the 

m = T~ = T-, = T„ case. In particular, the percent contribution 
2,w 3,w 1 4 7 
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of the integrations along the characteristic lines bd, ab, and he and 
of the reflected wave at point c are indicated therein. (The reflected 
wave at point c represents the contribution of all the characteristic 
lines in figure 2 other than lines hd, ah, and he, as is mentioned in 
appendix F . ) For weak shocks, the major contribution to shock attenua- 
tion comes from the integration along line hd. With increasing Mg, the 
contribution of the characteristic line he increases gradu ally to a val- 
ue of about 30 percent at Mg = 6.0. The integration along line ah leads 

to compression waves (which tend to accelerate the shock) and therefore 
is negative in figure 6 . Its value decreases to about -20 percent at 
Mg s 2.0 for the laminar case and to about -35 percent at Mg = 2.75 | 
for the turbulent case, and then increases with increases in M g . (The 
influence of characteristic line ab may be somewhat overestimated in 
the present analysis because of the assumption of an expansion me of 
zero thickness . ) The contribution of the reflected wave at point c var- 
ies from a value of zero at Mg = 1 to about -5 percent at Mg = 6. , 

The neglect of the latter contribution appears reasonable for the range 
of Mg considered herein. 


RESUITS AND DISCUSSION 

Equations (6) and (7) define the attenuation in a shock tube when 
the boundary layer is wholly laminar or turbulent, respectively. These 
equations are now compared with the experimental and theoretical results 
of reference 6. 


Comparison with Experiments of Reference 6 

Measurements of shock attenuation were obtained in the investiga- 
tion of reference 6 by using a high-pressure shock tube having a l/6- 
by l/8-foot rectangular cross section. Air at room temperature was used 
in regions 1 and 4. The air in region 1 was maintained at atmospheric 
pressure. Four sets of runs, corresponding to p = 4.061, 5.764, 

7.455, and 17.915 are reported. The results of these runs are summar- 
ized in figures 7(a) to (d), respectively. Some of the theoretical 
curves of reference 6 are included in these figures. The theoretical 
predictions of equations (6) and (7) are also indicated in figure 7. 

The latter were found from figure 3 with p-^ = 2117 and by assuming 

the wall to be a perfect conductor (since the shock tube had metal walls 
except for a pair of schlieren glass inserts). 

The Reynolds number at the contact surface for these tests is given 
directly by figure 3. The boundary -layer thickness at point b of region 
2 for the test conditions of reference 6 is plotted in figure 8. It is 
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assumed in figure 8 that the boundary layer in region 2 is wholly turbu- 
‘ lent and that the wall is a perfect conductor. The boundary -layer thick- 

ness is presented in the form of a relative thickness 25g b /h where h 
is the smallest shock-tube dimension normal to the flow and equals l/& 
for the shock tube of reference 6 . When 28g -b/h is small compared to 
1 , the assumption of a thin boundary layer is valid, and the theory of 
the present report is applicable. It may be seen from figure 8 that 
26g b /h < 0.20 for x ^ 12 for the test conditions of reference 6 . 
g Since 28g b /h tends to overestimate the effective boundary-layer thick- 

5 ness (because of the asymptotic manner in which the velocities approach 

free-stream conditions at the edge of the boundary layer), figure 8 in- 
dicates that the theory of the present report is applicable for compar- 
ison with the experiments of reference 6 . 

In figure 7(a), the first data point (x = 3') is close to the the- 
oretical value for a wholly laminar boundary layer. The other points 
fall somewhat above the theoretical curve for a wholly turbulent bound- 
ary layer. The Reynolds number per foot at the contact surface for this 

case is Reg b /x = 0.4x10® (fig. '3). Using Re = 0.5KL0® as a rough 

indication of the transition Reynolds number, as discussed in appendix 
H, it is reasonable to have the first data point near the laminar curve 
and the other points near the turbulent curve. 

The experimental results for Mg = 1.442 are given in figure 7(b). 
The Reynolds number at the contact surface is Reg t)/x = 0.8x10®. The 
data agree quite well with the values for a wholly turbulent boundary 
layer. There is a -slight tendency for the points to lie above the the- 
oretical curve which might be attributed to the short length of laminar 
boundary layer directly behind the shock wave. 

The experimental results for Mg = 1.518 (Reg - 5 /x = l.l x 10®) sire 
given in figure 7(c) and are in excellent agreement with the turbulent 
boundary-layer theory of the present report. 

In figure 7(d), the experimental results (Mg =* 1.792, 

Reg b /x = 2 . 8 > 10 ® ) fall somewhat below the theoretical predictions of 
the^ turbulent -boundary-layer theory of the present report. 

In general, the theory seems to agree reasonably well with experi- 
ment for the range of data considered in figure 7. Figure 7(d) indi- 
cates the poorest correlation and suggests that the present theory may 
underestimate the attenuation corresponding to large values of M g . 

properly computed displacement thickness would probably give a 
better estimate of the effective boundary-layer thickness. 
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Comparison with. Theory of Reference 6 

The present analysis assumes a relatively thin boundary layer, and 
the calculations are based on the vertical velocity at the edge of the 
boundary layer. The numerical value of v depends on terms which are 
related to the boundary-layer velocity profile, dissipation, and heat 
transfer. See, for example, equations (D3a) and (D3c) . These terms 
correspond to the use of wall shear, dissipation, and heat transfer in 
reference 6. Dissipation and heat transfer play similar roles (quali- 
tatively) in the present analysis and in reference 6. However, it can 
be shown that the velocity -profile term in the equation for v has a 
sign opposite to that of the wa ll -shear term in reference 6 for charap- 
teristic lines ab and be (e.g., appendix I) . Therefore, the integra- 
tions along lines ab and be, in reference 6, overestimate and under- 
estimate, respectively, the contributions of these characteristic lines 
to shock attenuation. If the same boundary -layer theory is used, the 
attenuation calculations of the present report and of reference 6 should 
be in quantitative agreement only in the limiting case Mg -*• 1 (for 

which the contributions of characteristic lines ab and be are negligible) 
and for those values of M. where the errors due to characteristic 

o 

lines ab and be tend to compensate. 

Since different boundary-layer theories were used in reference 6 
and herein, it is not directly possible to separate discrepancies be- 
tween the two methods due to the respective boundary-layer theories from 
the discrepancies due to basing the attenuation calculations on wall- 
shear, dissipation, and heat-transfer terms rather than on v. An 
estimate of the latter discrepancy can be obtained by reversing the sign 
of the first term in the equations for Lg (eq. (D3c)) and (eq. 

(D4c)) for the integrations of v along characteristic lines be and ab. 
Such a procedure shows that the theory of the present report and that 
of reference 6 agree at Mg = 1. With increasing Mg, reference 6 first 

overestimates the attenuation because of the increasing importance of 
characteristic line ab. At Mg = 0(l.5), reference 6 overestimates at- 
tenuation by about 10 to 15 percent. With further increases in Mg, the 

errors in characteristic line be become important and tend to compensate 
for the errors in line ab so that the discrepancy between the two methods 
decreases. At Mg = 0(2), the two methods are again in approximate 

agreement. As M increases further, line ab becomes relatively 
s 

less important compared with be, and the theory of reference 6 underes- 
timates the attenuation. At Mg = 6, the method of reference 6 appears 

to underestimate the attenuation by about 50 percent. The above fig- 
ures are only appro ximat e, because of the manner in which they were ob- 
tained, but they indicate the proper trend with Mg. 
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Prom figure 5(t>) it can be seen that the attenuation theory of ref- 
erence 6 (for turbulent boundary layers) is in good agreement with that 
of the present report for 1-^ Mg< 2. .The agreement is somewhat better 
than that to be expected from the discussion of the previous paragraph 
and is due, in part, to the use of an incompressible -boundary-layer 
theory in reference 6 . Hence, either theory could be used to correlate 
the experimental data of figure 7 . For Mg > 2, reference 6 consider- 
ably underestimates the attenuation. The laminar-boundary -layer theory 
of reference 6 differs by a factor of about 2 from that of the present 
lo report; therefore, there is a large discrepancy between the theoretieal 

° laminar -boundary-layer curves In figure 7(c). 


COHCLUDING- REMARKS 

A method is presented for computing the attenuation of a shock wave 
due to unsteady-boundary-layer action. The various assumptions involved 
in the analysis are summarized in this section, since these define the 
limitations of the method and suggest possible fields for improvements: 

1. Small perturbations . The equations of motion were linearized 
assuming the potential flow external to the wall boundary layer under- 
goes only small perturbations. For long shock tubes with large amounts 
of attenuation, it might be advisable to employ a characteristic method. 

2. Thin boundary layer. The assumption of a thin boundary layer 
(relative to shock-tube diameter) is consistent with the assumption of 
small perturbations of the potential flow. If the viscous effects span 
the entire tube cross section (i.e., long tubes), it may be advisable to 
base the shock -attenuat ion theoiy on wall shear. 

3. One -dimensional longitudinal waves. It is previously pointed 
out that the perturbations in a shock tube should be computed by assum- 
ing each element of wall surface area to be an elemental acoustic source 
of strength proportional to v. This gives rise to a complex wave pat- 
tern involving both longitudinal and transverse waves. For the purposes 
of the present analysis it was assumed that the longitudinal waves are 
of primary interest and that these can be computed on the basis of a 
simplified one-dimensional theory. This assumption is accurate when an 
observer is relatively far frcm the sources (since the details of the 
source distribution around a perimeter then become relatively unimpor- 
tant) but introduces errors when the observer is near the sources. It 
would seem that the assumption of one -dimensional longitudinal waves is 
accurate for weak shocks but introduces errors for stronger shocks, par- 
ticularly when the flow relative to the wal 1 is supersonic. 
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4. Boundary -layer theory . The boundary-layer theory of reference 

2 -was used. It may be assumed that the laminax-boundary-layer solution 
is reliable except possibly for very strong shock waves . In the latter 
case, very large temperature gradients exist normal to the wall, and it 
may be advisable to choose a different reference temperature from that 
used herein. Also, for strong shock waves, dissociation might occur 
which would also require changes in the laminar -boundary-layer theory. 

The turbulent-boundary -layer theory of reference 2 requires experimental 
verification for even the weak-shock case. However, the good agreement 
between the attenuation calculation based on the turbulent boundary layer 
of reference 2 and the experiments of reference 6 suggests that the 
turbulent-boundary-layer theory of reference 2 gives reasonable results, 
at least for the weak-shock case. 

5. Expansion wave of zero thickness. In order to simplify the prob- 
lem of determining the boundary layer behind an expansion wave, it was 
assumed in reference 2 that the expansion wave is of negligible thick- 
ness (i.e., "expansion shock"). This assumption is valid for weak ex- 
pansion waves but is in error for strong expansion waves. The boundary- 
layer solution for region 3 might be Improved for the strong-wave case 
by considering the finite thickness of the expansion wave. However, the 
contribution to shock attenuation of region 3 becomes small for the 
strong-wave case, so that an improved boundary -layer solution for region 

3 may not significantly affect the attenuation calculation. 

In addition to the preceding discussion, the following extensions 
of the present report might be pursued: 

1. The details for obtaining the perturbations at a fixed point in 
the shock tube, as opposed to finding the shock attenuation as is done 
herein, might be treated. This would be useful for further correlating 
the theory with experiment and for determining conditions at an aero- 
dynamic model (when the shock tube is used as a wind tunnel) . 

2. Formulas equivalent to equations (6) and (7) might be derived 
for the case where the boundary layer is partly laminar and partly 
turbulent . 

3. The bound ary-layer theory as outlined in appendixes D and E 
might be simplified, particularly for the turbulent case, so as to give 
reasonably accurate results without too tedious a boundary-layer 
calculation . 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 24, 1956 
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APPENDIX A 


SYMBOLS 

The following symbols are used in this report: 

A cross-sectional area of tube 

a speed of sound 

C eq. (C3) 

Cp specific heat at constant pressure 

D eq. (C3) 

d hydraulic diameter, 4A/l 

E eq. (C7 ) 

P eq. (C7) 

h smallest shock-tube dimension normal to flow 

appendixes D and E 

perimeter of shock-tube cross section 
Mach number of flow relative to wall 
shock Mach number relative to wall 


h’ L 5 

l 

M 

Ms 

MgjMs 

m 

V n 3 

P 

Re 

Re 2,b 


Mach number of flow in regions 2 and 3, relative to wall 

/' 

rate of mass addition per unit area per unit x 

appendixes D and E 

pressure, lb/sq ft 

Reynolds number (appendix H) 

Reynolds number at point b of region 2, 

2 


aix 




appendix H 
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T temperature, °R 

T r temperature of insulated wall 

t time 


u velocity of flow relative to wall 

u velocity of shock wave relative to wall 

s 

u^u^ velocity in regions 2 and 3 relative to wall 

v vertical velocity (positive when directed into tube) at 

edge of boundary layer 

x longitudinal, distance 

T ratio of specific heats 

A perturbation quantity (dp = perturbation of p, etc.) 

&2 laminar -boundary-layer thickness in region 2 (eq. (D5)) 

turbulent-boundary-layer thickness in region 2 (eq. (E3)) 
p coefficient of viscosity 

v kinematic viscosity 

£ integration variable representing x 

p . mass density 


a Prandtl number 

integration variable representing t 
wall shear stress, positive in +x-directicn 

Subscripts: 


l, 2,3,4 
a,b,c,d,e 

m. 

w 


regions of shock tube (fig. l) 

points on characteristic lines (fig. 2) 

evaluated at mean temperature of boundary layer 
(appendixes D and E) 

evaluated at wall 


4 
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Superscripts: 

+ associated with wave moving downstream (+x-direction) 

associated with wave moving upstream (-x-direction) 
Special notation: 

a 32 = a 3 /a 2 j p 12 = Pp/P2j etc. (Two successive integer subscripts, 
not separated by a comma, represent a ratio.) 


to 

» 

u 

o 
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APPEMDIX B 


CHARACTERISTIC -LINE GEOMETRY AND SOME INTEGRALS 

The equations of the characteristic lines considered in the present 
report are summarized. Seme important integrals are indicated. It is 
assumed that the characteristic lines are straight (corresponding to the 
ideal-flow case) and that the expansion fan can he considered as an 
"expansion shock" (i.e., expansion wave of zero thickness) moving with 
velocity a^ into region 4 (following sketch) . The point d with coor- 
dinates x,t is assumed to he on the shock-wave characteristic. 
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The equations of the characteristic lines are: 
Line hd 

g = x - + u 2)(t - t ) 




Line ah 


Line he 


t = t - 


X - 4 


a 2 + u 2 


u s T - £ 


|~1 + M 2 - (u, 

L 1 + M 2 


s/ a 2)l 


(x - €) 


^ - ( a 3 + u 3^K - 


T = T 


gh - g 
h " a 3 + u 3 




a.T + ? = 


1 + M 3 + a 43 


1 + 


“3 


(g -g a ) 


€ = ?h + ( a 2 " u 2)C' t b ” T ) 


U S T - ? = 


€- ^h 

T = T- h 

h a 2 - u 2 


Ti - M 2 + {ug/a^)! 

L 1 -“Z J 


fe e -g) 


(Bl) 
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In terms of x,t, tbe coordinates of points a, b, and c are 


Point a 


*a = - X 


Point b 


<c b 


Point c 


' C c = t 


Some important line integrals are 

.x 

r 

1 


[l+Mg - 

(u s / a 2 )^ 

a 4 

L i + m 3 

+ a 43 J 

u s 

[l + Mg - 

(V a 2)l 


L i + m 3 

+ a 43 J 


[i + m 2 - 

( u s/ a 2)J 

u 2 

u s 

[l + Mg - 

KA 2 )] 


fl + Mg - 

C u s/ a 2 ) 

-i 

[l - M 2 4 

KAg). 

- 

fl + Mg - 

( u 8 M 

- 

b- - Mg 4 

hsM. 





J 


~ n 2 


l-n c 


14^ I K T ~ %) (l-u 2 /u s ) (l + ^ 

^b 


(■ 


(B2) 




r\%C 

hi 


J 


b 

£b 


^ - «’ n2 - 1 


-n 3 


-n?f l~n 2 \ 
x 

A " n 2 j 


n 2 f 1 " n 2 


R) R-S)V 


n 2 > 


(a 4 T+ C) 34 = (l+U3/a 4 ) 


i-nj -s|)4 

1 + M 3 + a 43 \1 - n 3 y 




(B3) 
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APPENDIX C 


DJTERACTION OF PRESSURE WAVES WITH INTERFACES 

Weak pressure waves are assumed to overtake a contact surface or a 
shock wave . The strength of the reflected and transmitted waves is 
found. 


Contact Surface 

Consider the contact surface which separates region 3 from region 
2 (following sketch) . At a certain instant, known incident waves Ap t 

o • D 



and intersect the surface. The problem is to find the final 

waves Ap“ ^ and Apg -5 • From isentropic flow relations, 

"X. 

= pa Au* 


Ap" = -pa Au.' 

The boundary conditions across the contact surface are 

+ 


M,b + ^,b = ^J,b + ^3,b 
*4,b + ^2,b = ^3,b + ^i,b. 


(Cl) 


(C2) 
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From equations (Cl) and (C2), the expressions for AgvJ ^ and Ap^ 
are 3 3 


where 


^2,b “ “^b + “®s,b 
^3,b ~ r 32 a 23 n ^2,b ~ V± ^Z,\i 


C ~ ^ r 32 a 23 " 1 ^/CT 32 a 23 + x ) 
D ~ 2 /^ r 32 a 23 + 


(C3) 


The quantities C and D are reflection and transmission coefficients, 
respectively. 


Shock Wave 

Consider a shock wave moving with velocity u g into a stationary 
fluid (region l) . The region behind the shock is designated as region 
2. Conditions in regions 1 and 2 can be found from, shock-wave theory. 
In particular, with = x^} 

2ri»4 - Cri - I 

p 2i ~ n + 1 
- = -4 t(k 

a i r x + 1 \ B 

_L. 

At a certain instant, a known wave Ap 
ing sketch) . 3 




intersects the shock (follow- 
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» 

I 


The problem is to determine the net perturbation of the shock -wave. 
From equations (C4), 


AP2,a 

*2 

a l 


4r x 

T-l + 1 p 12 M s^ M s,d 





(C5) 


The net pressure and 
hy 


velocity perturbations in region 2 are then related 

m| 

^2,3 = 2p l a l *>2,4 CC6) 


With 3» 2 ,d = fip 2,d + Ap 2,d and Au 2,d ~ *4^1 + *4,d< us±n B equations 
(Cl) gives 

(C7) 

where 


^,a “ ^z.a 

^2,3 = FA 4 ,a 


E - ^12*12 x ^m 2 " l ( 2Pl2Sl2 1 !"m 2 + *) 

F s ^>i2 a i2 rr^j 77^ + 


The quantity E is a reflection coefficient. The resulting perturba- 
tion of the shock Mach number is 


^,d r l + 1 P 21 „ ^2,d 
M s “ M^ p 2 


(C8) 
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A ?2 d 

For Mg = 1, — = 0- For Mg infinitely large. 


Ap' 


^2,d V ^1 


2,d 

1 


APjjj d 

2‘ 


fri - i 


+ l 


Apg d 

, which becomes *— = - 0.14 for T-, 

^2,d 


- 1.4, 


indicating that the reflected wave is relatively small and is opposite 
in sign to the incident wave. 


if 


t?r*r\T, 
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APPENDIX D 


LAMINAR BOUNDARY LAYER BEHIND WAVE 


The laminar boundary layer behind a shock or thin expansion wave 
is analyzed in reference 2. Some of the results are s umma rized herein 
It is assumed that the wall temperature behind the wave is constant and 
that a and c^ are independent of temperature . The fluid properties 

jjl and k are referenced to a mean temperature as discussed in appendix 
C of reference 2. The notation of . reference 2 is compared with the no- 
tation of the present report in the following equations. The left-hand 
side of each equation represents the notation of reference 2, while the 
right-hand side is in the notation of the present report. 


Region 2 


Sr = 


= u. 


u e = U_ - 


s 

s "2 


Sr = 1 

u e 1 " (-VV 


A 


> 


v e ’ v 2 


x = UgTT - £ 


Sr , S>/ U s 

Ug ■ 1 - K/Ug) 


(Dl) 


J 


Region 3 


Sr 


= a 


u = a, + 
e 


'4 

4 + u 3 


Sr = 1 • 

u e 1 + (ug/aj 




v e = v 3 


x = a 4 T + £ 


Sr 

— - 1 = 


-Uj/a^ 


> (D2) 


u„ 


1 + (u3/a 4 ) ) 


Vertical Velocity at Edge of Boundary Layer 

The expression for the vertical velocity at the edge of the bound- 
ary layer can now be written as follows: 


Region 2 


. r 4k f 'S y 2 t 

2 ^ [_1 - lV u s)J \ u s T - 


(D3a) 
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\A lere 

n 2 = l/2 ( D3b ) 



-0.50+0. 47 (u 2 /u g ) 




(D5h) 

(D3i) 


The above quantities can he evaluated if the wall surface temperature 
T 2 w is known. The value of T 2 ^ depends on the heat capacity and 

conduction properties of the wall' and 'the wall thickness. A method for 
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evaluating is presented in reference 2. From the discussion in 

reference 2, it may he concluded that in most cases T 2 w approximately 
equals when the wall is a metal. ? 


Region 3 


where 


u|/ a 4 V n3 f v 3 ^ 

’> ■ % rife 




(D4a) 


n 3 = 1/2 


(D4h) 


II ^5,m 

|2 n 3 .i. 


T 3 f f(Tl-f)(l + u 3 /a 4 )“| 


T 3-1 


% 


•43 


* € 


r 3 + 


1 + ( u 3 / a 4 ) f T 3,w . T 5, 

V a 4 ' V T 2 



IT S 3 *1 


0 


(D4c) 


lim 


("H-f ) Cl + u 3 /a 4 ) 

U3/a 4 


1.217 


1 + (u^a^) 
2.326 + (u^/a^) 


'0 


0 • 22 u 3 / a 4 

1 + (u^a^) l+(u 3 /a 4 ) 

r 3 d,, = 1.686 /|j 5 - m - rT — ^ C„ 3 ) 


(D4d) 


(D4e) 


£ 


-0 . 50-0 . 36 (u-j/ a 4 ) 

1 + (u^a/) 1+ (u^a^) 

■ 1 -* 17 ti g.giS + (VH> ( ° S> (Mf> 


r 3 (°) = (a 3 ) 
T 


0.37+0. 50^5/ a 4 ) 

l+(u3/ a 4 ) 


3,1 ' T3 ^M|r 3 ( 0 ) 


T 


= 1 + 


(D4g) 

(D4h) 


2 
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ilf = °' 5 (^f + $ + °' 22 (^ ' x ) (M1) 

Equations (D4) can "be evaluated if T 3 w is known. If the wall is in- 
itially at temperature T^ then T 3 w = T 4 is generally a good esti- 
mate for T 3 w , as is discussed in reference 2 . 


Boundary -Layer Thi ckness in Region 2 


In the body of the report it is assumed that the boundary layer is 
thin relative to the shock-tube diameter. It is therefore of interest 
to present an expression defining boundary-layer thickness . Region 2 is 
of particular interest since the boundary layer in region 2 is generally 
larger than that in region 3. Let represent the laminar -boirndary- 

layer thickness in region 2, defined so as to correspond to 
From equation (18) of reference Z, 


— = 0.99. 
u 2 


5 2 = 



^2,m T 2 B 

v z( u b x ~ 

2.89V1 - i^/Ug T 2 - 1 w2 r 

2 , 

i-SSl 

u a < 

^ T 2,m 1 

U s 

LVI.413 - VS 2 h J 0 Z 



(D5) 


(Eq. (D5) uses a mean reference temperature and an interpolation formula 
for Hg (defined in ref. 3) and takes the upper limit on the integrals 

to be 00 rather than T]g.) To evaluate S 2 along the characteristic 


( 


1 + M2 - u s/ a 2 


1 + M 2 


line bd of figure 2, take u g T - £ equal to 
The value of 5g at point b is found by replacing u g T 

li 




U - K) 


K 


fl + M 2 - u s/ a 2 V, u 2\ 

\ 1 + “2 A u sJ X ' 
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APPENDIX E 


TURBULENT BOUNDARY LAYER BEHIND WAVE 

The turbulent boundary layer behind a shock or thin expansion wave 
is also studied in reference 2 for the case of constant wall surface 
temperature . The correspondence between the notation of reference 2 
and that of the present report was previously noted in equations (Dl) 
and (D2). The turbulent -boundary-layer solution of reference 2 (relat- 
ing to v) is now summarized. It is assumed that the wall surface tem- 
perature is essentially constant . 


Vertical Velocity at Edge of Boundary 'Layer 


Region 2 


i l-n«- 


v 2 = 


where 


L 2 = 0.0460 


n 2 = 1/5 
&2 [(u2/u s ) - l] 

5 VV u s ) 


$2 


(v - «) 

(Ela) 

- 

(Elb) 

. S 2 (V u s) 
bEVV - X ]J 

(Elc) 

3 ’*■■) - (‘ - 3)1 

(Eld) 


Q _2 W^s" 1 ) 

8 2 K/ U s) 


7 ’ 2 

T 

(l + Mi + "2, 

T 2,w 

•*■2,6 " 

V \J 2 > 7 % ^ 8 J 


1 - CV u s) 


1/4 , n 3/4 

2 




(Ele) 


(Elf) 
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'0 


X 


T- 


h* 

To 


£ dz 

2 

(Elg) 

( 1 + b 2 z - c 2 z 


T 2 r ' 

T 2,v 

(Elh) 


(Eli) 

-VW'* 

(Elj) 

! + 1 ) +0 - 22 (-if - 1 ) 

(Elk) 

essentially equals T-^ 

can be estab- 


^2,m 

To 


lished by the methods of reference 2. The assumption T 2 ^ w - T^_ appears 
reasonable for most cases. 


The integral Ig N is tabulated in reference 2 for N = 6, 1, and 
8 and various values of b 2 and c 2 . The integral can be evaluated 
analytically from the following expressions: 


where 



(Ell) 


(Elm) 

(Elh) 


(Elo) 
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^,11^ " X 2‘ 


N 


In 


f V ~ 1N \ i /jJ 

\ ^2 ) + l—/ m 1^-2 

hp=1 


m 


Alternate expressions for (I2 j^) and (ig 8X6 


fi ^ = 1 V f 1 Y 

^ 2,W^2 1- P 2 Z-y (tol) : \1 - P 2 J 

m=0 w 

" L 1 + U^ui-PaJ + ^^7 


+ • 


(I 2,N^ 2 (K+l)(l-\,) 


(N+lJCl-^) 


1 - (N + 1) (N+m.) (N-nm-1) (0 

m=l V 

, (W+l) (fl+l) 

[ ". Cn+1)(K+2)X 2 " ( w + 2 )( ih . 3 )X 2 - 


Region 3 


v 3 = L 3 


r u i/ a 4 


— il-n. 


1 + U s /a 4 ) 


(- 2 S -) 1 

\ a 4 T + V 


where 


n_ = 1/5 


L 3 = 0.0460 


S 5 |j- + ^ u 5/ a 4 f £3(113^4) \ ^ 


^3^ U 3/ a 4^ \ 0 3^ + ( u 3/ a 4 f| 


8 3 L 1 + ^3/ a 40 

6 3 (V* 4 > 


H 


u* T 

1 + a - 7 T 
a 4 1 3 


T 3 f T ^3 _ \ 

*3,w V 3 ' 6 + a 4 5 ’ 7 ) 


(Elp) 

* (Elq) 

1 ^ 

> 

_ J 

(EIx) 

(E2a) 

(E2b) 

(E2c) 


(E2d) 
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8 5 [l 4- (u3/a 4 )] _ 7( T3 /T 3|g ) V ( V \ “ 3 T 1 , 

8 5 ( U3 /a 4 ) ' 1 + tVa 4 ) |V ' V 1 ' a 4 J V ' ^ ^sj < E2e > 


/•„ ftA , „ yV * 

--(??) (ft) 


z K dz 


3,n / , , . 2 

70 1 + 3 3 2 " C 3 Z 


b 3 - C T 3,rA3,v) - i 
- [( T 3,r/ T 3> - *] <W> 


(E2f) 

(E2g) 

(E2h) 

(E2i) 


T 


= i + ! 


- 1 


•4 (°3> 


1/3 


T. 


5,m 

T, 




/t 3 

0.22 ° 


V 


T„ 


- 1 ) 


(E2j) 

(E2k) 


The conditions under which Tg w essentially equals T^ can he estab- 
lished by the methods of reference 2. The assumption T^ w = T^ appears 
reasonable for most cases. The integral 1^ ^ can be evaluated from 

the tables of reference 2 or equations (ElZ) to (Elr) (with the sub- 
script 2 replaced by the subscript 3 therein) . 


Boundary -Layer Thickness in Region 2 

The symbol 8g represents the turbulent -boundary-layer thickness 
as obtained frcm an integral (Karman-Pohlhausen) type of boundary-layer 
solution. From reference 2 } it can be shown that 


cc>r\T, 
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&2 



The value of 8 2 along characteristic line hd or at point b may be 
obtained by the substitutions indicated after equation (D5). 
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APPENDIX F 


DEVELOPMENT OF ATTENUATION FORMULAS 


The attenuation of the shock wave in a shock tube is now derived. 
The flow is considered to consist of the ideal shock-tube flow plus 
small perturbations due to the boundary layer. The ideal flow is uni- 
form in regions 1, 2 , 3, and 4, and is denoted by the appropriate sub- 
scripts (i.e., p^, p 2 , p^, etc.). Perturbations at any point are de- 
noted by A and an additional subscript indicating the point. Thus, 
the perturbations at point b of figure 2 evaluated on the right-hand 
side of the entropy discontinuity are designated by Ap 2 Aug and 

so forth. The expansion wave of the ideal flow is assumed to have neg- 
ligible thickness and to propagate into region 4 with the velocity 
u = - a^ as indicated in figure 2. (The same assumption was used in 
the attenuation study of ref. 6.) 

Let point d of figure 2 represent an arbitrary point on the shock 
wave. The problem is to find the net pressure perturbation behind the 
shock (i.e., Ap 2 . This requires an integration of equation (4) along 

r! 1 characteristic lines which contribute to Ap 2 Because of the 

entropy discontinuities, there are an infinite number of line segments 
along which the integration must be conducted (fig. 2). However, the 
major contributions to Ap 2 ^ can be shown to come from segments ab, 
cb, and bd. 


If Apt' represents the incident pressure wave at point d, then 

u • OL 


(J « u. 

equation (C7J gives 


But, from equation (4), 


^2,a = p ^2,a 




(FI) 


nX 


^2,d “ ^2,b + a 2 d(l . + .Mg) / V 2^^ T ^ 


(F2) 


'«b 


where the integration is conducted along the line bd. From equation (C3), 

(F3) 


Ap+ , = CAp^ , + IAp 


' 2,b 


'2,b 


.+ 

3,b 
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Again, from equation (4), 


2t2p 2 r c 

^2,b = A »i,c + a,a(l - Mg) / 

J «b 




> 


+ + ^ 3^3 

^3,b = ^3,a + a 3 d(l + Mg) 


v,(?/r)d£ 


(F4) 


'«a ^ 


It can be shown that Apt =0 in the present case. Moreover, 

O ■ 8. 

p 2 = p 3 . The expression for Apg ^ can then be expressed as 


1 ^2^ Ap 2,d 

F 2f 2 p 2 




In appendixes D and E it is shown that for wholly laminar or wholly 
turbulent boundary layers 



where L and n are independent of t and The values of L and 

n depend on whether the boundary layer behind a given wave is laminar 
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or turbulent . Substituting equations (F5) into equation (5) and using 
equations (B3) yield 



The corresponding change in shock Mach number can be found from equation 
(C5) . Appropriate values for L and n are given in appendixes D 
and E. 

If n 2 = n 3 and CE Ap 2 c /^9?2 d is assumed not to var y wibb x, 
then equation (F6) indicates that Ap 2 ^ c/^2, d = ^ and 



(F7) 


Equation (F7) indicates CE Ap g C /Ap 2 d to be independent of x and 

is therefore consistent, with the original assumption to this effect. 
Substitution of equation (F7) into equation (F6) gives accurate results 
even when Ug / n 3 , since the major contribution to Ap 2 ^ d comes from 

the first term on the right-hand side of equation (F6). The term 
CE Ap 2 C /Ap 2 d represents the contribution to Ap 2 ^ d of all the char- 
acteristic lines of figure 2 other than lines ab, be, and bd, and is re- 
ferred to as the contribution of the reflected wave at point c in fig- 
ure 6. It can be seen from figure 6 that this term is s ma l l and can 
probably be neglected in most cases . 
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If it is necessary to consider the "boundary layer "behind the shock 
(or expansion) wave as partly laminar and partly turbulent, then equa- 
tion (5) must he integrated accordingly. 
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APPENDIX G 

IDEAL SHOCK-TUBE RELATIONS 

- . Shock-tube relations, assuming ideal, flow, are presented herein 
for convenience. The formulas were obtained from reference 1. The 
notation p 12 = P-j/Pg, & 52 = and so forth > ls a S alQ used. 

Define 


a e (r + i)/(r - l) 
F s (r - i)/2y 


Then 


Pl4 p 21 


1 - (P 21 - 1) 




P4 ( c V^)p 

a l p 21 + 1 ^v ^*4 


p 34 “ ^ P 14 P 21^ 


1 hi 


1 + OqP2]_ 
21 4* P21 


Pot “ 


2 2P 4 

T 34 = a 34 = ( p 14 p 2l) 


T 21 “ a 21 ” 


h 


K + p 21 )~ 
+ a l p 21 J 


l/P. 


(Gla) 

(Gib) 

(Glc) 

(Gld) 

(Gle) 


U2 

°1 


u. 


P 21 


- 1 


P lW P 21 + !)] 1 2 


a 4 


1 - ( p 14 P 21> 




(Gif) 


(Gig) 
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^ s ^ = [pjU + c^)] / 

(Gill) 

u 2 _ 1L p 2l ■ 1 

^ ’ r l[jl p 2lW + p 2lj] 1/2 

(Gli) 

u 3 1 T "Pa 

-3 " ^ = I 4 r 4 L Cp i^ p! - X J 

(Glj) 
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APPENDIX H 


REYNOLDS NUMBERS AND TRANSITION 

A Reynolds number characterizing the houndary-layer development in 
regions 2 and 3 of the shock tube is now developed. 


Region 2 (— « -L < l\ 

\ u s u s T J ' 

In the notation of reference 2, a Reynolds number characterizing 



x ( u e -u w) 

the boundary layer behind a shock is defined therein as Re = 

v w u e 

Using the transformations indicated by equations (Dl) and arbitrarily 
basing v on free-stream conditions, the Reynolds number for region 2 
may be written as 


u 2 ( u s ^ - K) u^Ug 

Re 2 " v 2 (1 - W 


(HI) 


where are the coordinates of a point in region 2. Consider 

to be a point on the characteristic line through x,t (preceding sketch) . 
The Reynolds n umb er for points along this line, as a function of ^ is 
(using eq. (B2)) 


_ u 2 (x - g) f u 2 /u g VI + Mg - u s /aA 
Re 2 “ v 2 ~ «z/ u bA 1 + “2 / 


(H2) 


If the transition Reynolds number is known, equation (H2) can be used to 
determine the values of K at which transition occurs. The max imum 
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CO 

I 

O 

o 


Reynolds number in region 2 occurs at the contact surface (£ = and 
equals 



Equation (H3) can he used as a general index as to 'whether the boundary 
layer (along the characteristic line) is primarily laminar or primarily 
turbulent . A plot of Reg against Mg for air in region 1 (with 

T-^ = 520° R) is given in figure 6. 

The transition Reynolds number for the boundary layer behind a 
shock wave has not yet been established. In the absence of more accur- 
ate information, the transition Reynolds number for incompressible flow 

over a semi-inf inite flat plate ^Re = 0(0.5KL0 6 )j might he used to es- 
timate the transition point behind a shock wave. That is, take 
Reg = 0(0.5KL0®) as a rough estimate for the transition Reynolds number. 
This is probably a conservative estimate for the strong-shock cases, 
since the large amount of heat transfer to the shock-tube walls may have 
a very stabilizing effect on the boundary layer. 


Region 3 ( -1 ^ ^ — ) 

V a 4 T a 4 ) 

For a point in region 3, the Reynolds number as defined in refer- 

Region 



enee 2 becomes (using eqs. (D2) and v 


Re 3 - 


u 3 ( a 4 T + £) 


based on the free stream) 

V a 4 

1 + u 3/ a 4 


(H4) 
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For a point on the characteristic line influencing point x,t (preceding 
sketch ) , the Reynolds number becomes 

_ ^3^ u 5/ a 4 1 t Mg + a 45 f~j / 1 + M 2 - u g /a 2 \ a^j 

S 3 Vg 1 + Ug/a 4 1 + Mg jjc \^1 + Mg + a 43 ) uj 

At the contact surface (£ = %h)j equation (H5) has the value 

S %3 - *13 (sf) ( X + «2 - sf) M s . C®) 

The criterion Reg = 0(0 .5x10® ) might be taken as a rough estimate for 
the transition Reynolds number behind the expansion wave. 

The boundary-layer characteristics presented in appendixes D and E 
indicate a discontinuity in the boundary-layer profile across the contact 
surface. (Note that Re 2 -^/Reg ^ = Vg/v 2 .) This discontinuity does not 

actually occur, and the theoretical discontinuity thus represents a de- 
ficiency of the present method. The effect on shock attenuation is prob- 
ably small. 
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APPENDIX I 


GENERATION OF PRESSURE WAVES BY WALL SHEAR AND HEAT ADDITION 

The generation of pressure waves by body forces and. by heat sources 
in one-dimensional flow is investigated, and the results sire used to 
compare the method of the present report with that of reference 6. 


Wave Generation by Body Forces and Heat Sources 


Assume a uniform flow of pressure p, velocity u, and so forth, in 
a tube of constant cross section. Ibis uniform flow is assumed to be 
slightly perturbed by weak body forces and heat sources. The equations 
of motion are (neglecting the possibility of mass sources, since this 
case is treated in eq. (l)) 


P 


f SAu 


+ 



= f 


sap 

Sx 


cAp , SAii . cAp 

sr + p sr + u sr = 0 


> 


(n) 


cAs 

St 


+ u. 


cAs 

Sx 



J 


where f(x.t) is the body force per unit volume acting in the +x -direction 
and q(x,t) is the heat addition per unit volume per unit time. The sym- 
bol As represents the entropy perturbation of a particle and is related 
to the pressure and density perturbations by 


As _ Ap 
c v " P 



( 12 ) 


The perturbation at any point x,t can be shown to equal 


+ ^SL 
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The upper limit on the integral for Ap _ is +°> or -«> depending on 
whether M < 1 or M > 1. Since equations (l) and (ll) are linear, the 
solution for the case where mass sources are also present can he obtained 
by adding equation (2) to equation (13). Note that an elemental heat or 
mass source generates a sy mme tric pressure-wave pattern (positive pres- 
sure waves propagating in the downstream and upstream directions), while 
an elemental body force generates an antisymmetric pressure-wave pattern 
(positive pressure waves propagating in the downstream direction and neg- 
ative pressure waves propagating in the upstream direction) . 


Comparison with Reference 6 

In effect, reference 6 uses equations (13) to find the shock atten- 
uation in a shock tube. The value of q is obtained by averaging, 
across the tube cross section, the heat transfer at the wall and the 
viscous dissipation in the boundary layer. The value of f is found by 
averaging the wall shear across the tube cross section. Thus, if q w 
is the heat transferred into the boundary layer per unit wall area, 
is the net dissipation in the boundary layer per unit wall area, and T 
is the shear per unit area exerted by the wall on the fluid (taken to 
be positive in the +x-direction), then 


1 - 1 + V - 5 + V 


l - 4- 

f = A = d T w 


(14) 


Substituting equation ( 14) into equation (l3) and integrating along 
appropriate characteristic lines should yield the same results for shock 
attenuation as those which were obtained in reference 6 from a somewhat 
different viewpoint. However, these results are not in agreement with 
the results obtained from equation (4) (the latter being the basis of 
the present report). 
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Equation (4) is based on v while equation (l4) contains terms re- 
lating to the heat transfer at the wall, viscous dissipation, and wall 
shear. But, from, boundary-layer theory, it can be shown that v is 
dependent on terms related to the heat transfer at the wall, viscous 
dissipation, and the velocity profile (see, e.g., eqs. (D3a) and (D3c)). 
The heat-transfer and dissipation terms in equations (14) and (4) play 
the same qualitative roles, and therefore, trill not be discussed further. 
However, the wall-shear term in equation (14) in some cases has a sign 
opposite to that of the velocity-profile term in equation (4) . Thus, 
these terms are not always in qualitative agreement, and this leads to 
discrepancies between the results obtained from equations (14) and (4). 

Consider, for example, the boundary layer in a shock tube for the 
case of negligible heat transfer and dissipation. For this case, v de- 
pends only on the velocity -profile term. Similarly, only the x term 

is retained in equation (4). The signs of v and x^ and of the re- 
sulting pressure perturbations in regions 2 and 3 of the shock tube are 
summarized in the following table: 


Region 

Sign 


Perturbat ions 

Perturbations 


based on v 
(eq (4)) 

based on t. t 
( eq (14)) W 


V 

Ap + 

Ap~ 

X 

w 

Ap + 

£p~ 

2 

- 

- 

- 

■ - 

- 

+ 

3 

+ 

+ 

+ 

- 

- 

+ 


Thus, when heat transfer and dissipation are neglected, a perturbation 
solution based on v differs from a perturbation solution based on T w 
in regard to the sign of the Ap~ waves in region 2 and the Ap + waves 
in region 3. 

It may be concluded that attenuation solutions based on equation 
(14) differ from those based on equation (4) in the following respects: 
(l) The integration of equation (l4) along characteristic l i ne ab tends 
to overest ima te its contribution to shock attenuation, and (2) the inte- 
gration of equations (14) along characteristic line be tends to under- 
estimate its contribution to shock attenuation. For the limiting case 
-*■ 1, the integrations along ab and be are negligible, and equations 

D ' 

( 4 ) and (14) give the same results for shock attenuation. 
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When the boundary layer is thin, it is obvious that the attenuation 
calculation should be based on equation (4). If the viscous shear af- 
fects the entire cross section (as for a long shock tube), there is no 
longer a core of potential flow and the solution should probably be based 
.on equations (14) (or on the characteristics method of ref. 6 when non- 
linearities become important) . The quantities q^, and should 

then be based on unsteady pipe flow rather than on thin -boundary-layer 
theory. In practice, the boundary layer is generally sufficiently thin 
to permit the use of equation (4) . 
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Figure 3. - Reynolds number at point b for air in region 1 
at 520° R. Regjtj » v^gRe^ ^ (appendix H). 
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Norlral chock Mach number, 
(a) LmMnnr case (eq. (6)). 


Figure 6. - Attenuation, in air-air shock tube. >=> “ 520“ Prandtl number a , 0.70j ratio of specific heats T, 1.4. 
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Shock attenuation., (-1) 
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(h) Turbulent caaa (eq. (7)) 


Figure 5. - Concluded. Attenuation in air-air ebook tuba. T, ■= T. - 520° Ri Rrandtl nuabar a, 0.70 1 ratio 
of epeclfic heats t, 1.4. 
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(a) T.a.Tn1mr case. 


Figure 6. - Percent contribution to shock attenuation of 

characteristic lines of figure 2. Air-air shock tube; 

Tt = Oh = T„ = T_ = 520° R; Prandtl number a, 0.70; 
1 4 2,v 3,v ’ ’ 

ratio of specific heats x> 1.4. 
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(b) Turbulent case. 


Figure 6. - Concluded. Percent contribution to shock attenua- 
tion of characteristic lines of figure 2. Air-air shock 
tube; T x = T 4 = T 2jV = T 3jW = 520° R; Prandtl number a, 

0.70; ratio of specific heats X, 1.4. 
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O Experimental points (ref. 6) 

Theory, turbulent boundary layer 

(eg. (7) with 1 2fV = T S)V = T X ) 

Theory, laminar boundary layer 

(eq. (6) with T2,v ■= T3 ;V D T l) 
Theory, turbulent boundary layer 

(ref. 6, Cj o 0.0581 He”- 1 -/®) 
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Distance from diaphragm, x, ft 


(a) Pressure ratio p. /p , 4.061j shocfc Mach number relative to wall 
Mg, 1.344. 1 

Figure 7. - Pressure behind shock wave as function of distance from dia 
phragm. Air-air shock tube. T-j_ = T^ ° 520° R; hydraulic diameter 
d, 1/7. 
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(b) Pressure ratio p /p } 5.764; shock Mach number relative to wall 
M s , 1.442. ■ 


Figure 7. - Continued. Pressure behind shock wave as function of dis- 
tance from diaphragm. Air-air shock tube. Ti = T. = 520° R; hydraulic 
diameter d, l/7. * 


58 


MCA TIT 3278 


CM 


CM 

P< 


Q) 

> 


X 
O 
O 

XI 
tt> 

P 

•H 

X3 

0) 

rQ 


tn 

cn 

0) 

£ 



Distance from diaphragm, x, ft 


(c) Pressure ratio p /p , 7 ..455 ; shock Mach number relative to wall 
Mg, 1.518. 4 x 


Figure 7. - Continued. Pressure behind shock wave as function of dis- 
tance from diaphragm. Air-air shock tube. T n = T A = 520° R: hydraulic 
diameter d, l/7. 
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2 4 6.8 

Distance from diaphragm, x, ft 

(d) Pressure ratio p./p-,, 17.915; shock Mach number relative to wall 
M e , 1.792. 

Figure 7. - Concluded. Pressure behind shock wave as function of dis- 
tance from diaphragm. Air-air shock tube. Tj = T^ = 520° R; hydraulic 
diameter d, l/7. - 
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Figure 8- - Relative boundary- layer thickness at point b of region 2 for experimen- 


tal conditions of reference 6 (air in region 1; turbulent boundary layer; 
= T 2 ^ w ° &20° Rj h, l/B ft). 
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